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[57] ABSTRACT 

A method for formulating an exposure dose for an electron 
beam on a resist film for a pattern of geometric shapes which 
compensates for electron scattering effects utilizing hierar- 
chial design data which is preserved to as great as an extent 
as possible in the computation of the exposure dose. The 
exposure dose is corrected for both the forward scatter and 
backscatter effects of the electron beam in which the design 
data is modified for interactions of shapes which are affected 
only over the forward scatter range. In another version of the 
method, a multiple Gaussian approximation is used where 
the short term Gaussian terms are treated as the forward 
scatter terms and the long term Gaussian terms are treated as 
the back scatter terms.. 

9 Claims, 2 Drawing Sheets 
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ELECTRON BEAM PROXIMITY 
CORRECTION METHOD FOR 
HIERARCHICAL DESIGN DATA 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention generally relates to the technology 
of electron beam lithography, and more specifically to a 
method for proximity effect correction in electron beam 
lithography systems. 

2. Description of the Prior Art 

In recent years, various electron beam lithography sys- 
tems have been employed to write fine patterns on samples 
such as semiconductor wafers. The samples are typically 
covered by a thin film called a resist. The resist properties 
are modified by exposure to an electron beam according to 
a desired pattern. The patterned resist film is used as a stencil 
for subsequent processing. Electrons emitted towards the 
sample collide with nuclei, etc., of the sample and scatter. 
This is called backward scattering. A phenomenon called 
forward scattering also occurs. Forward scattering is the 
electron scattering that occurs in the resist prior to entry into 
the substrate. Forward scattering also represents the finite 
beam spread in real electron beam systems. This scattering 
of electrons may result in distortions to the desired pattern. 
This is called the proximity effect. The proximity effect may 
change the shape and dimensions of the part of the sample 
under the resist from those required. 

Various correction techniques to reduce the proximity 
effect have been developed. The most accurate techniques 
compute the effects of electron scattering on design patterns 
on a shape-to-shape basis, varying the electron dose applied 
to an individual shape, such that all shapes in a pattern 
receive the same total absorbed electron energy (i.e., the 
exposure). This is known as an exposure equalization tech- 
nique. 

An increase in the number of devices and elements on a 
chip increases the chip density and complexity. The chip 
density and complexity is further increased as features 
become smaller in size. This translates into greater demands 
on the computing resources required for performing prox- 
imity correction, as more shape-to-shape interactions must 
be computed. In addition, new electron beam lithography 
tools are now available which operate at high incident 
electron energies (approximately 100 Kev), increasing the 
range of back scattered electrons, resulting in an increase in 
the number of shape-to-shape interactions which must be 
computed. 

Common chip design tools employ a powerful data com- 
paction technique, known as nesting or design hierarchy, 
whereby a shape or a collection of shapes (i.e., cell) may be 
replicated any number of times in different parts of the 
design pattern without designing new shapes each time. This 
is done by simply referring to the original shape (or cell) for 
each subsequent placement, or instance, of that shape (or 
cell). 

In order to perform accurate proximity correction, it is 
necessary to know precisely the effects of all the shapes 
within the scattering range of a given shape. For a 100 KeV 
electron beam, this range may extend as far as 50 jum on a 
silicon substrate. Thus, the calculation of shape-to-shape 
interactions must be done for a large number of shapes. In 
general, this requires the loss of data compaction given by 
efficient use of design hierarchy, since the scattering effects 
influencing a given shape in a design cell may vary from one 
instance of the shape (or cell) to another. 
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In U.S. Pat. No. 5,051,598 issued Sep. 24, 1991, to Ashton 
et al., entitled, METHOD FOR CORRECTING PROXIM- 
ITY EFFECTS IN ELECTRON BEAM LITHOGRAPHY, 
incorporated herein by reference, a technique is described 

5 which provides a significant improvement in proximity 
correction efficiency, referred to as the dose-on-grid tech- 
nique. This technique is based on the splitting of the 
proximity function into two separate portions: a local 
portion, describing the effects due primarily to forward 

3Q scattering of the incident electron beam; and a remote 
portion, describing the effects due to long range back scatter 
of electrons. 

The back scatter portion is generally a more slowly 
varying function of distance from the incident electron beam 

15 than is the local portion, and it may be accurately repre- 
sented on a grid. According to the technique disclosed in 
Ashton et al., a fractional density matrix is generated. This 
matrix represents the pattern shape density on the grid, 
without reference to the individual shapes. A proximity 

2Q matrix is also computed, which represents the back scatter 
(remote) portion of the proximity function under a set of 
given electron beam exposure conditions (i.e., electron 
energy, substrate composition). When the density matrix is 
convolved with the proximity matrix, the resulting dose 

25 matrix, represents the back scatter corrected dose field on the 
grid. Individual shape doses are computed by interpolating 
the dose field from nearby grid points onto the shape. 

Despite the significant computational advantages offered 
by the dose-on-grid technique, the data compaction of the 

30 original design hierarchy is lost, since, in generating the 
fractional density matrix, the positions of all shapes in the 
pattern must be known. The result is twofold: (a) data 
storage volumes may increase dramatically as the fractional 
density matrix is generated by replacing individual shapes, 

3S and (b) in the calculation of the fractional density matrix, the 
same calculation may be repeated on the same shapes 
numerous times, resulting in inefficiency. 

For the foregoing reasons, there is a need for a method for 
efficiently computing a dose correction in electron beam 

40 lithography systems whereby the pattern design hierarchy is 
retained to as great an extent as possible during the process 
of proximity correction. 

SUMMARY OF THE INVENTION 

45 An object of the present invention is to provide a method 
for proximity effect correction in electron beam lithography, 
suitable for high beam voltages and/or for very dense 
patterns. 

Another object of the present invention is to provide a 
5 q method for proximity effect correction in electron beam 
lithography by dose correction. 

A further objective of the present invention is to provide 
a method of dose correction for reducing proximity effects 
in high voltage electron beam lithography wherein cells are 
55 imposed on a chip pattern and a correction process is applied 
based on fractional shape density matrices of patterns. 

Yet another objective of the present invention is to pro- 
vide a proximity effect correction technique for electron 
beam lithography including back scattering corrections and 
60 forward scattering corrections. 

Still another objective of the present invention is to 
provide a grid based dose correction technique for electron 
beam lithography which takes into account a proximity 
effect correction and a fractional density correction while 
65 retaining as much of the design hierarchy used in its calcu- 
lation as possible, thereby improving the computational 
efficiency of the calculation without sacrificing its accuracy. 
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Another objective of the present invention is to provide a the forward scatter range. These shapes or cells may either 

proximity correction technique which compensates for for- be retained as variants, or may have the dose information 

ward scattering without significant loss of the design hier- stored separately. For maximum efficiency, the assignment 

archy. of doses to individual shapes for both forward scatter and 

The present invention is directed to a method that satisfies 5 back scatter correction may be done in the electron beam 

these and other objectives. The method comprises, a com- pattern generation hardware. 

pact design hierarchy which is used to generate the frac- In another version of the present invention, a multiple 

tional density matrix. This may be done by traversing the Gaussian approximation to the electron energy distribution 

hierarchy tree of the design data and computing the frac- is used for improved accuracy. In this version the short range 

tional density matrix for the entire pattern coverage of basic *o Gaussian terms are treated as the forward scatter term, and 

cells (pattern cells containing actual shapes). As all instances the long range Gaussian terms are treated as the back scatter 

of a basic cell contain identical pattern coverage, the frac- term. 

tional density matrix may be computed very efficiently. Additional objects and advantages of the present inven- 

A point of novelty, among others, of the present invention tion will be set forth in the description which follows, and 

is that shape-to-shape interactions are explicitly calculated 15 in part will be obvious from the description, or may be 

only for the forward scatter correction, which is extremely learned by practice of the present invention. The objects and 

short ranged. For the vast majority of designs encountered in advantages of the present invention may be realized and 

electron beam lithography, the number of cells with different obtained by means of the instrumentalities and combinations 

environments over the forward scatter range will be much particularly pointed out in the appended claims. 

smaller than the number of cells with different environments 20 DDTCC nncpnTimnw nr ™c nn aauiktoc 

, , „ tt. u i 4* *■ BRIEF DESCRIPTION OF THE DRAWINGS 
over the back scatter range. The back scatter correction is 

calculated by a dose field technique, which involves the The accompanying drawings, which are incorporated in 

convolution of the proximity function for back scatter with and constitute a part of the specification, illustrate presently 

the pattern density matrix. Therefore, calculation of the preferred embodiments of the present invention and, 

density matrix may be done without flattening the pattern 25 together with the general description given above and the 

data. detailed description of the preferred embodiments given 

At the point where the back scatter correction dose is below, serve to explain the principles of the present inven- 

applied to the individual shapes, different instances of the U0D m wmcn: 

same shape which require different back scatter correction 3Q FIG. 1 is a flow diagram of the present invention, 

must be distinguished. However, since the correction has FIG. 2 is a typical hierarchial pattern layout in which the 

already been done, this step can be deferred until the time at method of the present invention is employed to determine a 

which the pattern is written by the lithography tool. This proximity correction. 

way, it is not necessary to store very large flat files. ncTATr _ n nccrDmnnM tuc 

j it _ . 4 . - 4 , , . , . , % DETAILED DESCRIPTION OF THE 

In another variation of the present invention, a local 35 PREFERRED EMBODIMENT OF THE 

density matrix may be calculated for basic cells. This may be INVENTION 

generated directly from the shapes in the cell, or 

alternatively, forward scatter correction may be performed Referring now in detail to FIG. 1, in which a method of 

on the individual shapes within the cell, yielding a weighted the present invention is illustrated, and generally referred to 

density matrix. This matrix is replicated over the entire 4Q by reference numeral 100. The method formulates an expo- 

pattern area, based on the instantiation (or placement) of the sure dose for an electron beam on a resist film for a pattern 

cell in the pattern. This is done for all cells in the pattern, of geometric shapes compensating for electron scattering 

resulting in a fractional density matrix for the entire pattern. effects, in which the design data is provided in step 110 in 

The back scatter correction is computed as described by hierarchial form. The design data consists of a library of cell 

Ashton et al., which is incorporated herein by reference, by 45 l yP es specifying the geometric shapes and their respective 

convolving the fractional density matrix with the proximity cells. 

matrix to create a dose matrix. The back scatter corrections The hierarchy is then searched in step 120 to identify 

are then assigned to the individual shapes in the pattern. shapes from different cells which interact over the forward 

Here too, the correction can take advantage of the design scatter range, act. The interacting shapes are referred to as 

data hierarchy, since it is necessary to apply different dose 50 intrusions. The searching step 120 determines if intrusions 

assignments only to those instances of a given shape or cell exist in the hierarchial pattern layout, 

for which the dose field varies. The dose required to expose The design hierarchy is then modified in step 130 by 

individual shapes may be stored as part of the individual either moving the intruding shapes into the cells with which 

shape information, in which case, variant shapes (or cells) they interact, or by creating variant cells in which the 

must be generated for those shapes (or cells) which have 55 presence of the intruding shapes is recognized. For typical 

variations in the individual shape dose assignments. designs, this will be required for only a very small percent- 

Alternatively, for greater data compaction, the dose infor- age of the total shapes. 

mation may be stored separately, with variant doses for The step of modifying 130 the hierarchial design data 

shapes (or cells) having variations in the individual shape typically comprises shapes processing which includes union 

dose assignments, but the original design hierarchy may be 60 of shapes, removing overlaps between shapes, expanding or 

restored, if desired. This requires that a reliable technique be shrinking shapes, fracturing shapes (cutting large shapes 

provided to reassign the correct doses to the appropriate into smaller shapes for optimized processing), and/or apply - 

shapes at the time of electron beam exposure. ing sleeves to large shapes. This sort of shapes processing is 

Forward scatter correction may be done on a shape-to- standard in typical data preparation for electron beam lithog- 

shape basis within individual pattern cells, with variants 65 raphy. 

created for instances of shapes or cells having shapes which An example of modifying the hierarchy design data in 

interact with shapes in closely neighboring cells, i.e., within step 130 is shown in FIG. 2 in which a typical hierarchial 
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pattern layout is illustrated. Cells A, B, C, and D are all The procedure used in generating the fractional density 

different from one another. Al, A2, and A3 are instances of matrix in step 160 is done by traversing the hierarchy tree of 

cell A; Bl, B2, and B3 are instances of cell B; CI, C2, and the pattern library, and computing the fractional shape 

C3 are instances of cell C. coverage for each hierarchial cell in the design (i.e., the 

The placement of cells A, B, and C is such that they are 5 fraction of the total area of the cell covered by shapes). For 

within the back scatter interaction range, C p, of one another, C f\ wh f* caa f^ C n J*"* - m ^ .^T^' !u" 

but not within the forward scatter interaction range, aa, densi y 0 Qe f* ceils *» projected directly onto the 

* e ii ™ u* u • * * n n ,u parent cell. The fractional shape coverage for the prune cell 

except for cell C3, which interacts with cell D over the K. tU «. . 4 . , / m . .? .J c t . 

f , . (i.e., the highest level parent cell) is the projection of the 

forward scatter range. shapc coveragc of all underlying cells in the hierarchy tree. 

The shapes in cells Al, A2, and A3 may be corrected for 10 a mesh of square cells is then superimposed onto the area of 

forward scatter as native cell A; the shapes in cells Bl, B2, the prime cell. The area of each cell is p/b, as described by 

and B3 may be corrected for forward scatter as native cell B; Ashton et al., which is incorporated herein by its reference, 

the shapes in cells CI and C2 may be corrected for forward The fractional shape density for each cell in the matrix is 

scatter as native cell C. However, cell C3 must either be then computed. Optionally, the fractional density may be 

combined with cell D, or be treated as a distinct variant from 15 weighted by the forward scatter corrected doses applied to 

native cell C, which may be referred to as native cell C\ the shapes in each cell. 

Thus, the hierarchy is retained for cells A1-A3, B1-B3, The proximity matrix (P) is then computed in step 170 

and C1-C2. There is only a slight increase in the pattern data precisely as is described in Ashton et al., which is incorpo- 

for variant cell C\ r&ted herein by its reference. The proximity matrix gives the 

Referring back to FIG. 1, after the hierarchy is modified 20 interaction of cells based on their relative geometrical con- 
in step 130, proximity correction for forward scattering is figuration and the proximity function, which, in the present 
computed in step 140. The forward scatter dose correction is invention, is only the back scatter term. The proximity 
computed on a shape-to-shape basis. As pointed out by matnx & computed in step 170 using one of a number of 
Ashton et al, which is incorporated herein by reference, any 25 Hnear combination proximity correction algorithms known 
desired or convenient existing algorithm may be used. * n tne art * 

The critical difference in the present invention is not in the ^ dose matrix (P) is then computed, 180, by convolv- 

correction algorithms, but rather in how the shapes are lD g the pactional density matrix (F) and the proximity 

grouped for modification in step 130. Since the forward matnx ( P ) to c° me U P Wlth a dose for each cell, 

scatter range is extremely short (of the order of the minimum 30 A back scatter correction is then determined in step 190 

feature size), for the purpose of computing the correction for from tne results of the previous three steps (160, 170, and 

any individual shape, interactions with other shapes must be I 80 )- |* * s determined in step 190 for each shape by 

considered only over a range -aa, which is the range of averaging the dose matrix (or dose field) over the shape by 

influence of forward scattered electrons. * polynomial or suitable interpolation function. 

Here, a is the Gaussian parameter used to characterize the 35 Alternatively, the back scatter correction may be determined 
forward scatter portion of the electron energy density dis- ™ ste P 190 on m area weighted average basis, 
tribution in a resist film, and is a multiplicative factor, Lastly a final corrected dose is computed by either inter- 
typically -2-3. For most shapes in a typical pattern design, polating the back scatter doses onto the shapes and corn- 
there will be a small number of interacting shapes within this fining the results with the forward scatter doses in step 
range. Moreover, the overwhelming majority of these inter- 40 200fl > or alternatively by assigning the back scatter dose 
acting shapes will be in the same cell in the design hierarchy. correction to the pattern areas and adjusting the exposure for 
Shapes from nearby cells which fall within the forward th e dose matrix cells by the back scatter dose in step 200£>. 
scatter range are called intrusions. Minor modifications to From the foregoing, it becomes readily apparent to one 
the hierarchy may be made in step 130 to account for these skilled in the art that the novel Electron Beam Proximity 
intruding shapes, such as moving the shapes to another cell 45 Correction Method of the present invention offers advan- 
in the hierarchy, or creating a variant cell to include the tages over the prior art electron beam proximity correction 
intruding shapes. The forward scatter dose correction may methods. Due to the inventive retention of the pattern design 
then be computed in step 140 for the shapes in each of the hierarchy during the process of proximity correction, the 
cells plus the variants. For most pattern designs which advantages offered by the inventive method resides in: 
employ hierarchial compaction, this will result insignificant 50 (a) increased computational efficiency in that multiple 
savings relative to flattening all the cells and calculating the instances of shapes with identical forward scatter envi- 
correction for each individual shape instance. ronments but different back scatter environments can 

After computing the proximity correction for forward be corrected only once for forward scatter; 

scatter in step 140, the forward scatter dose is stored by (b) reduced data storage requirements relative to flat file 

either assigning it to the shapes as an attribute to be retained 55 handling; and 

in the hierarchical design in step 150a or, alternatively, ( c ) optimized accuracy achieved by employing different 

writing it to a separate forward scatter dose file in step 1506, algorithms for forward scatter correction and back 

which may be combined with the shapes at a later point. scatter corrections. 

The fractional density matrix (F) is then generated in step While there has been shown and described what is con - 

160 as described by Ashton et al., which is incorporated 60 sidered to be preferred embodiments of the invention, it will, 

herein by reference, by superimposing a mesh of square cells of course, be understood that various modifications and 

onto the pattern area. The size of each cell is p^ (Ashton et changes in form or detail could readily be made without 

al. uses (3/a to describe the matrix cell size). However, the departing from the spirit of the invention. It is therefore 

difference between the present invention and Ashton et al. is intended that the invention be not limited to the exact forms 

that in the present invention the fractional density matrix 65 described and illustrated, but should be constructed to cover 

may be generated in step 160 directly from the hierarchial all modifications that may fall within the scope of the 

pattern library, without flattening the data. appended claims. 
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Having thus described our invention, what we claim as 
new, and desire to secure by Letters Patent is: 

1. A method for formulating an exposure dose for an 
electron beam on a resist film for a pattern of geometric 
shapes which compensates for electron scattering effects, the 5 
method comprising the steps of: 

providing design data in hierarchial form, the design data 
consisting of a library of a plurality of cell types 
specifying the geometric shapes and their location 
within the respective cell, 

searching the design data to identify shapes or cells which 
interact over a forward scatter range, the interacting 
shapes or cells being intrusions, 

modifying the hierarchial design data to account for the 15 
intrusions, 

computing a forward scatter dose correction on a shape- 
to-shape basis for the modified hierarchial design data, 

storing the data indicative of said forward scatter dose 
compensation, 20 

generating a fractional density matrix from the hierarchial 
library data, where the matrix grid size is a function of 
a back scatter range, 

computing a proximity matrix based on the back scatter 
range, 

computing a dose matrix by convolving the fractional 

density matrix and the proximity matrix, 
determining the back scatter dose correction based on the 

fractional density matrix, proximity matrix and dose 30 

matrix, and 

computing the final corrected dose for each shape in the 
pattern, said final corrected dose utilized by an electron 
beam pattern generation tool for patterning geometric 
shapes in a resist film with reduced electron scattering 35 
effects. 

2. The method as claimed in claim 1, wherein the modi- 
fying step comprises one of removing overlaps between 
shapes, union of shapes, expanding shapes, shrinking 
shapes, fracturing shapes, and applying sleeves to shapes. 40 

3. The method as claimed in claim 1, wherein the step of 
storing the forward scatter dose correction comprises assign- 
ing the forward scatter dose correction to the shapes as 
attributes to be retained in the hierarchial design data. 



25 



4. The method as claimed in claim 1, wherein the step of 
storing the forward scatter dose correction comprises writing 
the forward scatter dose correction to a forward scatter dose 
file to be combined with the shapes. 

5. The method as claimed in claim 1, wherein the step of 
generating a fractional density matrix comprises the sub- 
steps: 

traversing the hierarchy design data for the pattern of 
shapes, 

computing a fractional shape coverage for each cell in the 
design data, 

projecting the density of nested cells in the hierarchy 
design data onto a parent cell, the highest level parent 
cell being a prime cell, 

computing an overall fractional shape coverage by pro- 
jecting the shape coverage of all underlying cells in the 
hierarchial design data, superimposing a mesh of 
square cells onto the area of the prime cell, and 

computing the fractional shape density for each cell in the 
matrix. 

6. The method as claimed in claim 5, further comprising 
the step of weighting the fractional density by the forward 
scatter corrected doses applied to the shapes in each cell. 

7. The method as claimed in claim 1, wherein the step of 
computing the final corrected dose for each shape in the 
pattern comprises: 

interpolating the back scatter correction onto the shapes, 
and 

combining the interpolated back scatter correction with 
the forward scatter dose. 

8. The method as claimed in claim 1, wherein the step of 
computing the final corrected dose for each shape in the 
pattern comprises: 

assigning the back scatter dose correction to the pattern 
areas, and 

adjusting the exposure for the dose matrix cells by the 
back scatter dose correction. 

9. The method as claimed in claim 1, wherein the forward 
scatter dose corrections are treated as short range Gaussian 
terms and the backscatter dose corrections are treated as 
long range Gaussian terms. 
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